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Single quantum emitters like atoms are well- 
known as non-classical light sources which can 
produce photons one by one at given times[lj, 
with reduced intensity noise. However, the light 
field emitted by a single atom can exhibit much 
richer dynamics. A prominent example is the 
predicted ability for a single atom to pro- 
duce quadrature-squeezed light [3|, with sub-shot- 
noise amplitude or phase fluctuations. It has 
long been foreseen, though, that such squeez- 
ing would be "at least an order of magnitude 
more difficult" to observe than the emission of 
single photons j^. Squeezed beams have been 
generated using macroscopic and mesoscopic me- 
dia down to a few tens of atoms 0], but despite 
experimental efforts 04^, single-atom squeezing 
has so far escaped observation. Here we gener- 
ate squeezed light with a single atom in a high- 
flnesse optical resonator. The strong coupling of 
the atom to the cavity field induces a genuine 
quantum mechanical nonlinearity[9], several or- 
ders of magnitu de larger than for usual macro- 
scopic media[lCI H12l |. T his produces observable 
quadrature squeezing [l3l 4l5t with an excitation 
beam containing on average only two photons per 
system lifetime. In sharp contrast to the emis- 
sion of single photons |16| . the squeezed light 
stems from the quantum co here nce of photon 
pairs emitted from the svstem|l7j. The ability of 
a single atom to induce strong coherent interac- 
tions between propagating photons opens up new 
perspectives for ph otonic quantum logic with sin- 
gle emitters (isMiat . 



Unlike in a standard Kerr medium, our squeezing does 
not result from a simple nonlinear polarization of the 
medium but from a cavity-enhanced atomic coherence 
which exists for weak coherent driving. Consider a two- 
state atom with ground and excited states \g) and |e). 
In the absence of a resonator, the amount of squeezing 
is governed by the atomic coherence, cr = |5)(e|, and 
the excited-state occupation probability. The latter pro- 
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duces incoherent scattering which destroys the squeez- 
ing. Therefore, the laser intensity must remain low to 
preserve the atomic and hence the optical coherence, see 
Supplementary Information. Under this condition, the 
optical squeezing is determined by the fluctuations of the 
atomic coherence, Acr^ = {{a — (a) )^), itself given by 
Acr^ = — ( (T )^ owing to the fermionic character of a two- 
state atom, = 0. Note that this sets an upper bound 
to the amount of squeezing which can be obtained, even 
when all the light scattered by the atom in all directions 
is observed. 

The presence of the cavity introduces two important 
ingredients as sketched in Fig. 1. First, the cavity mir- 
rors spatially direct the squeezed light towards the de- 
tectors thus eliminating the need to observe the full Att 
solid angle. Second, the strong coupling to the opti- 
cal cavity mode makes the energy-level structure of the 
atom-cavity system anharmonic and thus allows for two- 
photon transitions |17| between the system ground state 
and the second dressed-state manifold containing two en- 
ergy quantaQ, |2±), sec Fig. la. The net result is an 
amount of squeezing at the output mirror that is given 
by —K{a)'^, where K depends only on the frequencies 
and width of the second manifold |2± ), see Supplemen- 
tary Information. It follows that for a given excitation 
and coherence of the atom, the squeezing using a cavity 
is scaled by the factor K which can be large for strong 
coupling. In particular, for a resonant excitation of the 
dressed states |2± ), K increases with the coherent atom- 
photon coupling rate, g, relative to the total rate of deco- 
hcrcncc (2k for the two-photon coherence and 7 for the 
atomic coherence decay), \K\ ~ g/i^K + 7). Wc empha- 
size that squeezed light produced by a single atom in free 
space would be anti-bunched ji^ . In our case, the light 
is squeezed and bunched [l7j. 

Our transition scheme is similar to that of a four- 
wave mixing process [2^, but the underlying physics is 
radically different as the scheme arises from the strong 
coupling of the quantized cavity field with a two-level 
atom. Moreover, the non-linear process differs from that 
in microwave experiments [2^ where short unitary evolu- 
tions interrupted by measurements produce non-classical 
field states while in our case the squeezed light is gener- 
ated and propagated out of a dissipativc resonator under 
steady-state driving conditions. 

The nonlinearity appears at a single-atom and single- 
photon level where quantum fluctuations play a major 
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FIG. 1: Four-photon process leading to quadrature squeez- 
ing for a nearly-resonant excitation of a single atom strongly 
coupled to the mode of an optical cavity. The cavity photon 
number states |0 ), |1 ), |2 )... and the atomic ground and ex- 
cited states \g ) and |e ) combine to form an anharmonic ladder 
of dressed states |n±) sharing n excitations. The system is 
excited with laser light of frequency uj. Two laser photons are 
absorbed (red arrows) and reemitted down the ladder of states 
(green and blue arrows) to produce squeezing. The temporal 
dynamics of this four-photon process as well as the spectrum 
of squeezing is revealed in a homodyne detection scheme, see 
sketch in (b): A high-finesse cavity containing a single ®^Rb 
is excited with a weak coherent beam. The transmitted pho- 
ton flux is monitored with a single-photon counter to control 
the atom-light coupling. The fleld properties of the reflected 
beam, picked up by an optical circulator (polarizing beam 
splitter and A/4 waveplate), are measured with the balanced 
homodyne detector. The phase of the measured quadrature 
is controlled with a piezoelectric actuator. 



role [HI, [lH so that its understanding requires a full quan- 
tum treatment rather than a simplified linearized ap- 
proach. In the experiment, the quadrature operator of 
the light field, Xg, is measured outside the cavity us- 
ing a homodyne detection with a controllable phase 6. 
In this way we measure the time and normally ordered 
(symbol ::) autocorrelation of the quadrature fluctua- 
tions AXg ^ Xg - {Xg), 



AXe(T)AX(,(0) 



MK{<^?f{r)). (1) 



Here 5R denotes the real part and the function /(r) de- 
scribes the dynamics of the emission process: the pho- 
ton pairs cascading via the states |1±), detuned by 
Ai-t = UJ — uji± with respect to the frequency uj of the 
probe laser, create beatnotes decaying according to the 
linewidths 71-1- of these states. Defining their (complex) 
detunings as u)i± = Ai-t -|- «7i± yields 



/(r) = a+ exp(iaji+T) + a- exp(iwi_r) 



(2) 



where a± depend only on iui± and sum up to one (/(O) = 
1). Finally, the spectrum of squeezing is obtained by a 
Fourier transform. 

Single ^^Rb atoms arc held inside a high-finesse optical 
cavity by a red-detuned 785 nm dipole trap (Fig. lb). A 



cavity mode, nearly resonant with a closed atomic transi- 
tion at 780.24 nm, is excited with a Pin = 8.5 pW coher- 
ent beam, containing on average 2.0 photons per cavity 
decay time (see Supplementary Information). The effec- 
tive atom-cavity coupling g/27r = 12 MHz exceeds the 
atomic dipole and cavity field decay rates (respectively 
7/27r = 3 MHz and K/27r = 1.3 MHz), bringing the sys- 
tem in the strong coupling regime. The coupling strength 
is verified by monitoring the transmitted light intensity 
using a single-photon counter. The quadratures of the 
light field refiected from the cavity arc measured with a 
homodyne detector and sampled with a high-resolution 
fast digitizer. After trapping and probing each atom, 
an additional reference data sample is acquired with an 
empty cavity, providing an accurate measurement of the 
shot noise level, and the phase of the local oscillator is 
shifted by ±7r/2 to alternate between the X = Xq and 
P = X^ji quadrature measurements. For each quadra- 
ture, we acquire « 3 s of strong-coupling and sa 30 s of 
reference data. 

A time-resolved acquisition provides direct access to 
the quadrature autocorrelations, revealing the dynam- 
ics of the atom-cavity system. For each quadrature 
Xg — X OT P we calculate the time-domain autocor- 
relations of the homodyne signal acquired for strongly 
coupled atoms, and subtract the autocorrelations of the 
empty-cavity reference, which leaves a quantity propor- 
tional to (: AXg{T)AXg{0) :). The normalization factor, 
obtained by measuring the mean value of the excitation 
field, perfectly matches the value expected from the de- 
tector's parameters and includes its overall 55% efficiency 
without artificially compensating any other experimental 
imperfection such as the losses in the atom-cavity system 
itself. 

Figure 2 presents two homodyne autocorrelation mea- 
surements. The first data set. Fig. 2a, is measured with 
the probe tuned on the empty cavity resonance, lu = lUc, 
while detuned from the trapped atom by 27r x 8 MHz. 
The data present oscillations with a 9 MHz frequency 
and a 50 ns damping time, characteristic of a beatnote 
with the closest one-photon dressed state |1+). The 
contribution from the state |1— ) is negligible, such that 
/(r) oc exp(icji4_T). The measured autocorrelations are 
clearly phase-dependent and the antisymmetry between 
the X and P quadratures, which translates into an an- 
tisymmetry of their noise spectra is a first sign of non- 
classicality. The value at r = corresponds to the differ- 
ence in the integrated noise variance between the signal 
and the reference: the negative value for the X quadra- 
ture confirms the presence of squeezing with one atom. 
For comparison, a second set of measurements, shown 
in Fig. 2b, has been performed with the probe tuned 
close to the two-photon resonance, (w — uJc)/2tt = —12 
MHz, while detuned from the atom by 27r x 3 MHz. In 
this case the probe frequency is closer to the state |1— ) 
(Ai_/27r = 9 MHz) but the transition through the state 
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FIG. 2: Photon beat in the time domain. Autocorrelation 
functions for the X and P quadratures. In a, left column, 
the probe is tuned on the empty cavity resonance, which 
favours the transition |2-|-) 19^), with a charac- 

teristic beat via state In b, right column, the probe is 

close to the two-photon resonance |2— ) with an interference 
between the two possible paths |2— ) |1— ) — >■ \gO) and 
|2— ) — i- — j> \gO). The theoretical curves are obtained 
from the analytical model described in the text, taking into 
account the extraction efficiency of intracavity photons and 
the dynamical response of the homodyne detector. Statisti- 
cal errorbars correspond to the standard deviations obtained 
with ~ 30 000 data samples and ~ 300 000 reference samples. 



(Ai+/27r = —18 MHz) has a comparable contribu- 
tion. With our parameters the two transitions interfere 
destructively, making the signal more complex. More- 
over, the resonant excitation at the two-photon transi- 
tion decreases the phase coherence and therefore, due 
to a larger contribution of the incoherent emission, the 
asymmetry of the X and P autocorrelations is slightly 
altered. The theoretical fits obtained from the analytical 
model above, including the extraction efficiency of intra- 
cavity photons, the dynamical response of the homodyne 
detector, and an additional 1 MHz decoherence rate due 
to atomic motion are in very good agreement with the 
experimental data for both parameter regimes. 

The noise spectra of each quadrature, relative to the 
shot noise extracted from the reference data, are deter- 
mined by a direct Fourier transform of the time-domain 
homodyne signal and corrected only for the 55% homo- 
dyne detection efficiency. They are presented in Fig. 3. 
For uj = ujc (Fig. 3a) at the expected 9 MHz « |Ai+|/27r 
frequency, we observe 12 ± 2 mdB of squeezing on the X 
quadrature, with the same amount of antisqueezing on 
P. With a "perfect" cavity, single-ended and free from 
absorption losses, this value would be 5 times larger. For 
Lu ^ uja (Fig. 3b) the negative interference between the 
two decay paths results in a cross-over between squeez- 
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FIG. 3: Squeezing from one atom. Noise spectra for the X 
and P quadratures relative to the reference empty-cavity sig- 
nal, corrected for the homodyne detector's efficiency, for the 
measurement of Fig. 2a and Fig. 2b. Each curve is an average 
of ~ 30.000 data and ~ 300.000 reference spectra obtained by 
Fourier-transforming the homodyne signal over 100 /is inter- 
vals (see text for details). The error bars correspond to the 
standard deviation and account for the statistical uncertainty 
and the systematic uncertainty of the shot noise level. The 
theoretical curves are obtained from the analytical model in 
the text, transposed into the frequency domain. 



ing and antisqueezing for each quadrature: the transition 
through the nearest state |1— ) leads to squeezing on X 
and antisqueezing on P at low frequencies, whereas the 
transition through |1-|-) leads to the opposite behaviour 
around 13 MHz. The spectra are in remarkable agree- 
ment with the analytical model of cavity quantum elec- 
trodynamics for both parameter regimes. 

It is instructive to compare the phase-dependent non- 
linearity reported here for one atom to those obtained 
with macroscopic media. For the measurements with 
uj = uJc (Fig. 2a and 3a), where the largest squeezing 
is observed, the mean intracavity photon number is only 
n = 0.033. Therefore, only a small part of the imping- 
ing beam experiences losses due to atomic spontaneous 
emission, mirror coating absorption and cavity transmis- 
sion, and T] = 86% of the input power P,;„ = 8.5 pW 
is reflected towards the homodyne detector. For a light 
beam weakly squeezed by a third-order non-linear pro- 
cess, we can write an input-output relation expressing the 
squeezing in the frequency domain {AX'^) = —r]rPin/2 
relative to the shot-noise level of 1/4, where r deter- 
mines the nonlinear response of the system. In our 
case rr] = 1.6 x 10^ W~^ which exceeds by 7 orders of 
magnitude the Kerr nonlinearity of a standard single- 
mode optical fibre[l2[ with the same amount of losses, 
and by 4 orders of magnitude the x^^^ nonlinearities ob- 
tained by four-wave mixing in macroscopic atomic sys- 
tems with similar band widths [l( 



111. The ~ 10 mdB 



squeezing level is of course very small compared to the 
10 dB achieved in state-of-the-art parametric upconver- 
sion/downconversion experiments [27[, which remain the 
best way to generate squeezing as a resource but require 
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several watts of pumping power: decreasing the latter 
to a '-^ 10 pW level would bring the squeezing down 
to ~ 10~^ dB. Finally, compared to experiments using 
atomic beams 0, trapping a single atom allows us to 
keep the non-linearity constant, as required for most ap- 
plications. 

The squeezing achieved in our setup is limited by the 
losses in the mirror coatings and the use of a symmetric 
cavity with two output ports. Both limitations will be 
removed in the near future by means of an asymmetric 
cavity with a lower loss rate. Other types of resonators 
with smaller mode volumes like microtoroidsjl^ or fi" 
bre resonators [i^ could allow one to generate squeezed 
light on an atom chip. Using artificial atoms like quan- 
tum dots in microcavitiesjlSf would lead to larger and 
fixed atom-cavity couplings. Furthermore, recent the- 
oretical and experimental progress, e.g. (30| indicate 
that such experiments could soon be transposed to mi- 
crowave systems using Josephson junctions coupled to 
strip-line resonators U 21 1 , where even stronger non- 
linearities could be achieved. 
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OBSERVATION OF SQUEEZED LIGHT FROM ONE ATOM EXCITED WITH TWO PHOTONS : 

SUPPLEMENTARY INFORMATION 



THEORY 

In this section we develop the theoretical description 
of the squeezing process in our atom-cavity system and 
outline the derivation of the equations used to interpo- 
late the experimental data. To complement previous ap- 
proaches based on the pure-state formalism [isl 31 1, we 



make a direct connection between the squeezing of the 
atomic dipole and the squeezing of the cavity field, and 
point at a distinction between the contributions of single- 
photon and two-photon excitation processes. 



Quadrature squeezing of the cavity field 

We consider a single cavity mode defined by the an- 
nihilation and creation operators, a and respectively, 
and obeying the canonical commutation rules [a, a^] = 1. 
The quadratures of the intracavity field are defined by 



1 



(3) 



where 9 is an adjustable phase. Defining the fluctuation 
of a with respect to the steady state value ( a ) as Aa — 
a— (a), the quadrature variance (AX^) = {{Xg — {Xg))^) 
reads: 

(AX«2> = lme-'''{Aa^)) + {Aa^ Aa)) + J, (4) 
which, when normally ordered (symbol ::), reduces to 

(: AXl :) = i(5R(e-2^^(Aa2)) + (AatAa)) . (5) 

For a coherent or a vacuum state, one has (Aa^) = 
(Aa^Aa) = 0, thus the quadrature variance is equal 
to the shot noise, (AX|) = 1/4, i.e. (: AX^ :) = 0. 
The state of the field will be quadrature-squeezed if for 
some 9 the variance of Xg drops below the shot noise, 
(: AX^ :) < 0. 

The term (AatAa) is the incoherent part of the spec- 
trum and can never be negative. In order to observe 
squeezing in a nearly-resonant dissipative system such as 
ours, this term must remain as small as possible, which 
can be achieved at weak enough excitation intensities. In 
this coherent (phase-sensitive) limit it becomes negligible 
compared to the coherent term (Aa'^) and the quadrature 
squeezing reduces to 



Squeezing generated by one atom in an optical 
cavity 

To calculate (Aa^), we consider a two state atom, with 
a ground state \g) and an excited state |e), interacting 
with the mode of the cavity with photon number states 
|0),|1),|2).... The relevant physical parameters are: 

• The atom-cavity coupling (7, 

• The cavity field decay rate k, 

• The atomic dipole decay rate 7, 

• The frequency of the driving laser w, 

• The amplitude of the driving laser e, 

• The detuning between the driving beam and the 
cavity resonance frequency Ac = w — Wc, 

• The detuning between the driving beam and the 
atomic resonance frequency A^ = cj — w^. 

The evolution of the system's density matrix p is given 
by a master equation 

-t: ^ ^P = -t[-H^./c + Hp,p\ + uLaP + -jLaP ■ (7) 
at n 

In the rotating wave approximation and in the interaction 
picture with respect to the laser frequency, the Jaynes- 
Cummings hamiltonian is 

Hjc/fi = -Acotfl - AoctV + g(aV + acr^), (8) 

where <t = \g){e\ and = |c)(3| arc, respectively, the 
lowering and raising operators. The terms 

KLaP = K{2apa'^ — pa^a — a^ap) (9) 
'jLcP ~ 'y{2apa^ ~ pa"^ a — ap), (10) 

respectively describe the losses through the cavity mir- 
rors and spontaneous emission from the atom into the 
free space modes. The laser light incident on the input 
mirror results in a coherent excitation of the cavity mode: 



Hp/h = e(at + a) 



(11) 



(: AX| :) « -^{e-^''^{Aa^)) 



(6) 



The time evolution of the expectation value of a time- 
independent operator O is computed as d{0)/dt = 
Tr{OCp). Defining complex detunings as Wc = Ac + 
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in, LOa ~ Aa + 17, we have 



with the result in the body of the paper, 



— (a) = i(5c(a) - - e) 
^(a^) = 2i((I)c(a^) - .g(acr) - e(a)) 



(12) 
(13) 
(14) 



— (acr) = i((a;a + a;c)(a(T) + g{a^(Jz) - e(o-)) (15) 

where = |e)(e| — |5)(.<7| is the population inversion of 
the atom. 

We operate in the weak excitation regime, where the 
population of the excited atomic state |e) is negligible. In 
this case (ctz) « —1, {aaz) ~ —(1)1 {o^o'z) ~ —{a^), and 
the equations above provide the steady-state solutions 



a = 



UJnUJr 



£.9 



(Aa^) 
(AaAcr) 



(l^c('^a - 


f Wc) - g'^)(^a^c 

t^g^^a-^^c) 




(t^c(tia - 


fWc) - g'^){u!aU!c 




the fluctuations: 






-e'g' 




(l^c(t^a 


-f Wc) 


~g^Y 


(£^0(^0 


-f Wc) ~ g^){uJaU!c 


-g^Y 



(16) 
(17) 

(18) 

(19) 

(20) 
(21) 



Notice that, under the assumption of weak excitation, the 
field amplitude (a) and the atomic coherence (tr) scale 
as e, whereas {<Y) and (acr) scale as e^. More involved 
calculations [15, 31 1 would show that (Aa^Aa) scales as 
e^. This is nothing else than saying that the spectrum is 
coherent at weak excitation, (a^a) « |(a)P oc e^, thereby 
justifying the use of Eql^Hlto compute EqlHl 

We now define the complex detunings of the dressed 
states \n±) as 



This allows us to simplify 
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,2\ _ T.^/ /_\2n 



(Aa^) = K{-{cjY), 
where the constant K is simply 

25^ 



K 



(22) 

(23) 
(24) 

(25) 



(:AX|:)=-i5R(e-2^^A'(a)2). (26) 



Physical interpretation of the squeezing mechanism 

Equation [H] shows that the squeezing of the optical 
field in the cavity is directly related to the squeezing of 
the atomic coherence. Due to the two-level structure of 
the atom, the elementary property cr^ = directly gives 
{{(7 — {a)Y} = —{<^Y- In the dressed-state picture, 
(u) depends only on the normal modes |1±), see Eql^H 
thus highlighting the one excitation processes. The con- 
stant K, in contrast, connects the squeezing of the cavity 
field to the squeezing of the atomic coherence and demon- 
strates the importance of the two-photon states |2±), as 
obvious from its definition Eq l25l The observed squeez- 
ing is therefore fundamentally dependent on the anhar- 
monic structure of the atom-cavity system. The under- 
lying physical mechanism can be summarized as such: 
At low atomic excitation, the dynamics of the internal 
state of the atom are largely coherent, which due to the 
quantized nature of the atom (two-level structure) di- 
rectly translates into polarization squeezing, — (a)^, and 
couples back to the cavity mode to create quadrature 
squeezing, —K{aY- 



Time-domain evolution and squeezing spectrum 

The quantum regression theorem provides a set of lin- 
ear differential equations for the two-time correlations: 



— (Aa(T)Aa(O)) = i(wc(Aa(r)Aa(0)) 

UT 



— (Aa(T)Aa(O)) ^ i(5,(Aa(r)Aa(0)) 

CtT 



.g(Aa(r)Aa(0))), (27) 
:D,(Aa(r)Aa(0)) 
5(Aa(T)Aa(0))). (28) 



Solved with the initial conditions given by Eqs. [20l and 
[2n it yields for t > 

(: AXe(r)AXe(0) :) = -^^R{e-^'^' K {af f (r)) , (29) 
where the regression of the fluctuations is given by 



fir) 



UJ1+ - iOl- 

5i_ 



- Wl- 



exp(ia;i_T) (30) 
exp(ia3i_|_T). (31) 



The squeezing spectrum measured outside the cavity, 
normalized to the shot noise, can be obtained by a simple 
Fourier transform of the autocorrelations. Taking into 
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account the overall detection efficiency rj and the fact 
that the fluctuations of the field leaking outside the cav- 
ity can be related to those inside by (Aa^^j) = 2K(Aa^), 
it reads 

/>oo 

Sg{n) = 1+7]2k / dTCOs(nr)(: AXg{T)AXg{0) :) 
Jo 

(f}2„i;2_ -^2_^2^)}' (32) 

sum of two Lorentzians with frequencies and widths de- 
fined by the normal modes. 

Interpretation in the dressed-state picture 

The physical description can now be transposed into 
the dressed-state picture to determine the transition 
mechanism shown on Fig. 1 in the main paper body. 
Two pump photons make a near-resonant excitation of 
the second doublet of dressed states, |2±), which decay 
via the normal modes |1±). The time dependence of 
the measured signal, determined by /(r), corresponds to 
the beatnote of these photons with the local oscillator 
which has the same frequency ui as the probe. It is the 
sum of two terms oscillating at the detuning frequencies 
oj — uJi± of the normal modes and damping according to 
their linewidths. For each of the two possible transitions, 
the two emitted photons appear in opposite sidebands 
with respect to the local oscillator and, due to the coher- 
ence of the emission process, this results in quadrature 
squeezing at the corresponding frequency. 

EXPERIMENTAL SETUP 

The atom-cavity system is based on a 123 ^m - long 
Fabry-Perot optical resonator formed by two identical 
mirrors (transmission 2.8 ppm, absorption losses 4.0 ppm, 
radius of curvature 200 mm, diameter 7.75 mm), reaching 
a finesse of ss 470.000. A 785 nm laser beam, resonant 
with a TEMqo cavity mode, is used to stabilize the cavity 
length and to form a red-detuned dipole trap for single 
®^Rb atoms, injected into the cavity by an atomic foun- 
tain. A 780.24 nm probing beam excites another TEMqo 
cavity mode, nearly resonant with the (5^5*1/2, F = 
3, mF = 3 — 5^P3/2, F = 4, = 4) atomic tran- 
sition. The probing power impinging on the cavity is set 
to 8.5 pW, which corresponds to 2.0 photons per rele- 
vant temporal mode defined by the cavity decay time of 
60 ns. The effective atom-light coupling g/2TT = 12 MHz 
exceeds the atomic dipole and cavity field decay rates 
(respectively ^/2tt = 3 MHz and K,j2-K — 1.3 MHz) and 
brings the system into the strong-coupling regime. The 



quality of the coupling, reduced compared to the maxi- 
mal value gmaxl'2'^ = 16 MHz by the motion of the atom 
inside the trap, is deduced from the cavity transmission 
measured with a single photon counter: strongly coupled 
atoms drive the cavity off resonance and decrease the 
transmitted photon flux by a factor of 50 @. The light 
exiting the cavity through the input mirror is picked up 
by an optical circulator and reflected towards a balanced 
homodyne detector (Thorlabs PDB120A). The phase of 
the local oscillator, controlled by two mirrors mounted 
on piezoelectric actuators, is locked using an auxiliary 
772 nm frequency-stabilized beam far off-resonant with 
respect to any cavity mode and hence directly reflected 
by the input cavity mirror. The homodyne signal is digi- 
tized with a 14-bit resolution at a 100 MHz rate. Within 
the overall 40 MHz bandwidth of the homodyne detection 
system, limited by the anti-alias filter of the digitizer, the 
800 yuW local oscillator power yields a shot noise level 
> 9 dB above the electronic noise. 



EXPERIMENTAL SEQUENCE 

The experimental sequence is divided in three parts: 
preparation, measurement, and control. During the 
preparation time (~ 3 s), atoms are loaded into a 
magneto-optical trap located 25 cm below the cavity, op- 
tically cooled down to w 5 /.tK, and launched towards 
the cavity which they reach at nearly zero velocity. A 
sharp drop in the cavity transmission heralds the arrival 
of a well-coupled atom, trapped by raising the dipole 
potential from 0.2 to 0.9 mK. This triggers a 20 ms mea- 
surement sequence, which depends on the probing fre- 
quency. When the probe beam is resonant with the cav- 
ity (w = Wc), a cavity cooling mechanism leads to good 
atomic storage times and the system is probed continu- 
ously with a Pin = 8.5 pW incoming power. When the 
probe is tuned close to the atomic resonance (w ~ Wa) no 
cooling occurs: in this case the measurement sequence 
is divided into alternating 200 /^s probing and cooling 
intervals (the probe power and frequency being set to 
Pin = 8.5 pW, w « Wa for probing and Pin = 1.7 pW, 
uj = Uc for cooling). After each measurement, the last 
part of the sequence (« 1 s) is used to verify the powers of 
the trapping and cooling beams, as well as the power and 
the phase of the local oscillator, and to shift this phase 
by ±7r/2 to alternate between the X and P quadrature 
measurements. 

The extremely low signal level combined with short 
atomic storage times puts drastic requirements on the 
experimental stability and technical noise suppression. 
This is achieved by actively controlling all experimental 
parameters and by optimizing the data acquisition proce- 
dure. During the 20 ms measurement sequence, the atom 
remains strongly coupled for 1.6 ms on average before it 
leaves the cavity. The remaining data, acquired when 
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FIG. 4: "Raw" homodyne autocorrelations, comparing the 
signal (strongly coupled atom) and reference (uncoupled 
atom) data for a given parameter set {X quadrature, probe 
frequency u = ljc, for other parameters the curves are very 
similar). The visible oscillations correspond to the impulse 
response of the detector, i.e. to the autocorrelations of the 
shot noise measured with the finite detection bandwidth. 



the cavity is empty, is used as a reference. This nearly 
simultaneous signal and reference acquisition leads to an 
excellent cancelation of technical noises and slow ther- 
mal drifts. For each probing frequency we repeat the 
measurement sequence during ~ 100 h to acquire ~ 3s 
of strong-coupling and « 30 s of reference data for each 
quadrature. 

DATA ANALYSIS 

For each trapped atom, we divide the 20 ms data sam- 
ple into 200 fi s intervals and, using the measured cav- 
ity transmission, we select those where the atom was 
strongly coupled (transmission T < O.OATmax), and 
those where the cavity was empty (T > 0.7 Tmax)- For 
each interval, we calculate the autocorrelations with the 
maximal 10 ns time resolution allowed by the digitizer's 
speed, we average the calculated functions over all quali- 
fied intervals with a given parameter set, and we subtract 
the mean value at large time delays {t > I /xs) where sig- 
nals become uncorrelated. 

Figure S] shows that the obtained "raw" autocorrela- 
tions of the strong-coupling signal are practically super- 
imposed with those of the empty-cavity reference. The 
oscillations visible at this scale correspond only to the 
impulse response of the detector, i.e. to the autocorre- 
lations of the shot noise measured with a finite band- 
width. Physically relevant signals, « lO'' times smaller, 
become only accessible by calculating the difference of 
these curves, where these large oscillations cancel out, 
leaving a quantity proportional to (: AXe{T)AXg{0) :). 



In order to properly normalize it, we unlock the cavity 
and measure the offset on the X quadrature: all of the 
light is then reflected towards the homodyne detector, 
which corresponds to measuring a coherent state \a) with 
{X) = a = •\/2l) and (P) = 0. The corresponding nor- 
malization factor perfectly matches the value expected 
from the detector's parameters and includes its overall 
efficiency rjd = 0.55, without artificially compensating 
any other experimental imperfection such as the losses in 
the atom-cavity system itself. 

An accurate subtraction between the signal and the 
reference requires an extremely stable shot noise level, 
determined by the local oscillator's power. This power is 
actively stabilized during the preparation and the control 
phases of the measurement sequence, but the stabiliza- 
tion system is too slow to follow the fast switching be- 
tween the 200^8 probing and cooling intervals, and dur- 
ing the 20 ms measurement time the power is controlled 
by a sample-and-hold circuit. This circuit presents a 
small systematic drift, and the shot noise level changes 
by « lOdBm (i.e. sa 0.2%) during these 20 ms. Since 
the strong coupling data corresponds to the beginning, 
and the reference data to the end of the acquisition, this 
variation appears on the differential autocorrelations, al- 
though it is too small to be visible on Fig. H) To measure 
and compensate for this drift, for each parameter set we 
select the sequences where the atom was not trapped at 
all (triggering the acquisition but leaving the trap within 
the first 100 /is). In this case the variance of the ho- 
modyne signal corresponds only to the shot noise. By 
averaging over all such sequences, we obtain the evolu- 
tion of the shot noise during the 20 ms acquisition. We 
then determine the average time interval when the atoms 
stayed strongly coupled (first 1.6 ms) and the one when 
the cavity was empty (last 16 ms of the 20 ms acquisi- 
tion) and, for each case, calculate the average shot noise 
variance. Rescaling the reference data by their ratio ^ 
(typically 0.998) compensates for the drift and allows to 
completely cancel out the contribution of the shot noise 
on the differential autocorrelations: dominant on Fig. |4j 
it is absent from Fig. 2 in the main body of the paper. 

Although the agreement between the measured 
data and the theoretical autocorrelation function (: 
AXg{T)AXs{0) :) (Eq. ^ is satisfactory, it can be fur- 
ther improved by taking into account the percussional 
response of the detector -D(t), obtained by renormaliz- 
ing to 1 the integrated shot noise autocorrelations on Fig. 
m The theoretical curves on Fig. 2 in the main paper 
body correspond convolution of the theoretical autocor- 
relations (: AXg{T)AXg{0) :) with D{t). 

We calculate the noise spectra with a 10 kHz resolution 
by a direct Fourier transform of the time-domain homo- 
dyne signal divided into 100 /xs intervals. Figure [5] shows 
that the "raw" spectra of the signal and the reference 
arc, like the autocorrelations, completely superimposed, 
the reference spectrum being less noisy due to better 
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FIG. 5: "Raw" signal and reference homodyne noise spec- 
tra {X quadrature, uj = Uc, for other parameters the curves 
are very similar). The visible features correspond to the fre- 
quency dependence of the detection gain, defined by the anti- 
alias filter (steeply rolling off beyond 40 MHz) , and to a small 
(200 mdB) and narrow (< 10 kHz) residual technical noise 
peak at 6.3 MHz with an equal height on both curves. The 
DC offset peak (-38.9 dBV) is cut off. 



statistics. The only visible features are the smooth vari- 
ations of the shot noise level due to a non-flat response 
of the digitizer's anti-alias filter, and a small (200 mdB) 



and narrow (< 10 kHz) residual technical noise peak at 
6.3 MHz with an equal height on both curves. We rescale 
the reference spectrum by the same factor ( as the auto- 
correlations to compensate for the small shot noise drift, 
and calculate the difference between the two spectra in 
logarithmic scale. Like for the autocorrelations, this can- 
cels out the residual technical noise, and brings the noise 
level for frequencies above 30 MHz within 1 mdB from 
0: as can be expected from the parameters of the atom- 
cavity system, no spectral features appear at those fre- 
quencies, and the noise level of the signal is equal to 
the shot noise level of the reference. Since the narrow 
10 kHz resolution is only necessary for efhcient technical 
noise cancelation and is otherwise excessive for a sys- 
tem where relevant decay times remain below 0.1 ^s, we 
average out narrow-band spectral fluctuations by convo- 
lution with a 1 MHz Lorentzian filter, which amounts to 
assuming that signals separated by t > 0.16/^s are un- 
corrclated. Finally, we correct the obtained spectra by 
the independently measured 55% homodyne detection ef- 
ficiency, to obtain the curves presented on Fig. 3 in the 
paper body. We verified that the same spectra can be ob- 
tained by Fourier-transforming the differential autocor- 
relation functions and rescaling the result by the cavity 
decay rate and the overall detection efficiency including 
cavity losses. 



